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Higgs hunting with B decays* 

Ulrich Nierste'^ 
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B physics is sensitive to the efltects of Higgs bosons in the Minimal Supersymmetric Standard Model, if the 
parameter tan/3 is large. I briefly summarise the role of B ^ M^M" ~^ t'^Vt in the hunt for new Higgs 

effects and present new results on the decay B —> DtUt-: Using the analyticity properties of form factors one can 
predict the ratio R = B{B — > DrVr) /B{B Dlui), £ = e,jj,, with small hadronic uncertainties. In the Standard 



Model one finds R = 0.31 ± 0.02, 



= (0.71 ±0.09)% and B{B" D+t- 



(0.66 ±0, 



/o, if 

the vector form factor of the Heavy Flavor Averaging Group is used. B —> DtUt is competitive with B^ t^Vt 
in the search for effects of charged Higgs bosons. Especially sensitive to the latter is the differential distribution 
in the decay chain B DvrT~[-^ -k'Vt]. 
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1. Higgs effects in B physics 

Weakly-coupled extensions of the Standard 
Model (SM) typically possess a richer Higgs sec- 
tor than the latter. The easiest extension of the 
SM Higgs sector involves one additional Higgs 
doublet and is realised in the Minimal Supersym- 
metric Standard Model (MSSM). At tree-level the 
MSSM Higgs sector coincides with a Two-Higgs- 
doublet model (2HDM) of type II, in which down- 
type fermions receive their masses solely from one 
doublet, while up- type fermion masses exclusively 
stem from Yukawa interactions with the other 
Higgs doublet. An important parameter is the 
ratio tan/3 of the two vacuum expectation values. 
In the type-II 2HDM the bottom and top Yukawa 
couplings Ub and yt satisfy the relation 
Vb mi, 



yt 



mt 



■ tan/3. 



Values around tan/3 = 0(60) correspond to yb^yt 
unification, which occurs in grand-unified theo- 
ries (GUTs) with a minimal Yukawa sector. The 
idea of grand unification seems to call for low- 
energy supersymmetry, which stabilises the elec- 
troweak scale against radiative corrections from 
heavy GUT particles, improves the unification of 
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the gauge couplings and reconciles the predic- 
tion of the proton lifetime with its experimen- 
tal bounds. Probing the large-tan /3 region of 
the MSSM is therefore of great interest, since 
the question of Yukawa unification sheds light 
on the Yukawa sector of the underlying GUT 
theory. Yet large values of tan/3 are also inter- 
esting from purely phenomenological considera- 
tions: The tension between the measured anoma- 
lous magnetic moment of the muon, [T], and 
the Standard Model prediction [5] invites super- 
symmetry with tan /3 10, and larger values of 
tan /3 allow to saturate with heavier super- 
partners. Recent global fits of electroweak and 
B-physics observables to the constrained MSSM 
and the model with minimal gauge-mediated su- 
persymmetry breaking gave best fits for values of 
tan/3 = 54 and tan/3 = 55, respectively [3]. 

B physics is excellently suited to study large- 
tan (3 scenarios, because down-type Yukawa cou- 
pHngs grow with tan/3 and tan/3 = 0(50) corre- 
sponds to 2/f, ~ 1 [1] . Most dramatic effects can be 
expected in the leptonic decays Bq ^+^^ (with 
q — d OT s and £ — e, fi or t), which are not only 
loop-suppressed in the Standard Model but also 
suffer from an additional helicity suppression. In 
particular the 95% CL limit 

B{B, fi+^-) <5.8 • 10"^ 

« 18-SSM(i?,->^V-) 



1 



2 



from the CDF experiment [5] already cuts into the 
large-tan /3 region of the MSSM parameter space 
[3]. This is even true for the popular scenario of 
Minimal Flavour Violation (MFV) [5], in which 
the supersymmetric contribution to the 
IJi^ li' amplitude involves the same elements of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
as the SM amplitude. For large tan (3 one finds 



2. Charged-Higgs effects at large tan/3 



cx 



(1) 



where M^o is the mass of the CP-odd Higgs bo- 
son and ey is a loop function which depends on 
several MSSM parameters [3]. The pattern of 
Eq. ([T]) is in sharp contrast with the case of the 
naive type-II 2HDM, in which B{Bs is 
proportional to only four powers of tan /? ^ . In 
the LHCb experiment it will be possible to mea- 
sure B{Bs fi^ iJL~) for any value of tan/3. One 
can then test the MFV hypothesis by checking 
whether B{Bd ^^fi^)/B{Bs A^^M ) agrees 
with \Vtd/Vts\'^f%Jfl^, where /s, is the decay 
constant of the Bq meson. 

Due to the overall loop factor of ey the bound 
on tan^ derived from Eq. ^ depends on 

a plethora of other MSSM parameters. By con- 
trast, effects of the charged Higgs boson i7+ en- 
ter B decays at the tree-level. The information 
gained from charged-Higgs-mediated processes is 
therefore more directly related to the parameters 
of the MSSM Higgs sector, with smaller depen- 
dences on e.g. superpartner masses. effects 
are best studied in leptonic and semi-leptonic B 
decays, in which hadronic uncertainties are un- 
der sufficient control. We specify our discussion 
to the case of a r lepton in the final state, because 
the Yukawa couplings of the third fermion gener- 
ation are largest. The B factories have observed 
the decay B tv with [5] 

B{B^Tv) = (1.41 ± 0.43) X lO"'^, (2) 

which allows to place first useful constraints on 
ian (3 /Mfj+ [5]. In the following sections I will 
elaborate on another promising chargcd-Higgs 
hunting ground, the decay B — )■ Dtv, and com- 
pare this mode with B — > tUt- The presented 
work has been performed in collaboration with 
Stephanie Trine and Susanne Westhoff |TD] . 



The qbH^ coupling (with q = u or c) is given 



by 



g rrib tan (3 



"''^^ 2^2 A% 1-1- eh tan/? 

q{\+^^)bH+., (3) 

where the small Yukawa coupling yq is set to 
zero. The bottom quark mass mf, is defined in the 
same QCD renormalisation scheme as the current 
5(1-1-75)6. In the 2HDM of type-II the parameter 
£(, vanishes. In the MSSM with MFV is a loop 
factor; the typically dominant squark-gluino con- 
tribution to £(, is proportional to ^* /Mg, where 
/i is the Higgsino mass parameter and Mg is 
the gluino mass d]. A priori ef, could be com- 
plex, but experimental constraints from electric 
dipole moments severely constrain the phase of 
fi*/Mg and thereby of Cb- Since |eb|tan/3 can 
be of order 1, the charged-Higgs phenomenology 
does involve genuine supersymmetric parameters, 
yet with much less impact than in B^ fj,~^^~. 
In the MSSM with a generic flavour structure €b 
is different for q ~ u and q — c and may ob- 
tain a sizable phase. In the generic 2HDM Cf, is 
generated at tree-level and is typically complex. 
For an early extensive study of B ^ Diiyg and 
B Drur in the MFV-MSSM see Ref. pUj- 

The effective hamiltonian describing h qrv 
transitions mediated by W'^ or iJ"*" can be writ- 
ten as 

iJeff = {97^(1 ~ 75)&T741 - 75);.. 

— q [gs + gpib\ ot(i - 75)i^r| 



h.c. 



(4) 



In the MSSM the effective couplings gs and gp 
read 



gs=gp-- 



tan2/3 



M]j+ {l + eb tan/3)(H-er tan/3)' 



(5) 



Here Et is the analogue of eb for the VttH^ cou- 
pling. B^ T^v probes the coupling gp [5]: 



B{B+^T+u) cx |K6r/i|l-5P 



(6) 
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Figure 1. B{B^ t^v) vs. gp for real gp. 

The B meson decay constant Jb — 216 ±38 MeV 
[T3| is the dominant source of theoretical uncer- 
tainty in the extraction of |1 — .gp| from B{B^ 
T'^v). Eq. ^ is confronted with the experimental 
result of Eq. © in Fig. [TJ 

The semi-tauonic decay B Drvr instead 
probes the coupling gs, because B and D have 
the same parity. In the MFV-MSSM B Drvr 
and B'^ — )■ t^v probe the same parameters, 
because in Eqs. ^ and ([5]) is the same for 
b u and b ~* c transitions. Thus within 
the MFV-MSSM we can combine the information 
from both decay modes to constrain the com- 
bination of tan /?, Mh+ , efc and in Eq. ^ . 
If new physics is found, a comparison of gs ex- 
tracted from B DtVt with gp obtained from 
5+ — > r+j/ will probe physics beyond the MFV- 
MSSM. B Dtv compares to B+ t+v as 
follows: 

i) B{B DrVr) ~ 506(B+ ^ T+i^r)- 

ii) \Vcb\ entering B Dtv^ is much better 
known than | Vub \ ■ 

iii) The uncertainty in lattice calculations of 
f]j needed for B — > rvr is 30%. Instead 
B — > Dtv involves hadronic form factors. 
We will see below in Sect. [3] that the associ- 
ated theoretical uncertainty is smaller, if ex- 
perimental data on _B ^ Dtv^ with £ = e, fi 
are exploited. 

iv) The three-body decay B DtVt has 
decay distributions which discriminate be- 
tween iy+ and H'^ exchange. 

v) i? — > TV.J- is mildly helicity-suppressed. 
In i? ^ DtVt the (transverse) 



contribution is helicity suppressed 
in the kinematic region with slow 
D meson (P-wave suppression) |14j : 

W-Boson: Higgs: 

B,D 



V V V 



0=1 s=i <=i 



<ZZ1 S=0 1=> 



3. Charged-Higgs effects in i? — > DrVr 

B — > Ddui, £ = e,ii,T, involves the hadronic 
matrix elements of the vector current and of the 
scalar current, which are expressed in terms of 
two form factors 111 Vi and Si: 



{D{pD)\crb\B{pB)) -~ 

Vi{w)—= 
2y/mBmD 



u . n 

Pb +Pd 



+ S'i(w) (1 + w)y/rriBrfiD 

{D{pD)\cb{f,)\B{pB)) = 

Si{w) (1 + w)^mBmD 



friB - mo 



ruB - mo 



, (7) 



where mo and ms are the meson masses, q = 
Pb ^ is the momentum transfer and the run- 
ning quark masses mt, and rUc are evaluated at 
the renormalisation scale /i at which the scalar 
current cb is defined. The kinematic variable 



2mp)niB 



(8) 



is defined in such a way that the kinematic end- 
point q^ = {niB — rnoY corresponds to w = 1. 
Heavy quark symmetry implies [15] 

= Fi(l) - l + 0(l/m„a,). (9) 

For mi = the other kinematic endpoint is 
at = corresponding to Wmax — (rn^ + 
m|,)/(2TO£)WB) = 1.59. The absence of a pole 



at = in Eq. ^ implies 



Si{w„ 



Vl(Wmax)- 



(10) 



The next step towards a precision analysis of 
B — > DtVt exploits the analyticity properties 
of form factors (TB]: The location of poles and 



t There is a typo in the definition of Si in Ref . I10| . 
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Figure 2. Vector form factor Vi{w). Dots: 
BELLE data on B ^ Dlvi [T^ with statis- 
tical errors only. Dark gray band: form fac- 
tor from BELLE data with HQET constraint at 
w — \ (systematic errors dominate at large re- 
coil). Light gray band: form factor from HFAG, 
which includes older CLEO and ATLAS data [15]. 
Solid Hue: best fit of BELLE data to and aX ■ 
Dashed line: best fit of Ref. [17] (using different 
parameters). 

branch points in Vi and 5*1 can be inferred from 
the particle spectrum. The conformal mapping 



^/{mB + moY - + \J (jtib + moY - to 

and the elimination of subthreshold poles renders 
the form factors analytic in the new variable z 
in the entire kinematic region (see |16|10j for de- 
tails). That is, we can paramaterise Vi and 5*1 in 
terms of power series in z. With a proper choice 
of the free parameter to the kinematic range 
1 < w < 1.59 is mapped onto < |z| < 0.032. 
B Dtvi with £ = e,fi only involves Vi. We 
can use experimental data to verify that the first 
two coefficients and aY of the power series 
are sufficient to describe the normalisation and 
shape of Vi. Moreover, the dependence on aY is 
moderate. The result is shown in Fig. [2l Finally 
we need the scalar form factor S{w). The corre- 
sponding parameters Og and af are fixed through 
Eqs. ([U and pO)) . using the l/rric and as cor- 
rections from |19j . An alternative approach uses 
lattice data to fix the form factors near w = 1 



|20j . Note that our analysis faces much smaller 
hadronic uncertainties than the extraction of | Vet \ 
from B Dlvf. We first fix |Vc6|Vi(w) from ex- 
periment. Then Eq. pop determines the normal- 
isation of |Vch|5i(u') in terms of measured quan- 
tities. Keeping only Oq already reproduces S{w) 
to 90%, and neither hadronic uncertainties nor 
the error in \Vcb\ have entered the prediction of 
B — !■ DtVt yet. Hadronic uncertainties only en- 
ter when we increase the accuracy further and 
fix the slope parameter af by including Eq. ([9]). 
These uncertainities are suppressed by l/nic- In 
this step one also needs the experimental value 
of \Vcb\, whose error is around 2%. A remaining 
source of hadronic uncertainty is the next term af 
of the series in z. In the data for in Fig. [2] we 
found the influence of aY negligible^ and expect 
the same behaviour for af and 5*1. 

The first application of our analysis in llOj is a 
new prediction of branching fractions in the SM: 



B{B- D"T-Dr) 

B{B° D+T-Vr) 

~ B{B Dtvi) 



(0.71 ±0.09)% 
(0.66 ±0.08)% 

0.31 ±0.02 



(11) 



using the HFAG form factor Vx [IH]. Note the 
small uncertainty in R which is the relevant quan- 
tity for charged-Higgs hunting. This has to be 
compared with the 0(40%) error of iKfjP/s en- 
tering the SM prediction of i? — > tv^-. The un- 
certainties in Eq. (|lip will decrease further with 
better data on i? — > Dii/g. The dependence of 
R on gs is shown in Fig. [31 A home-use for- 
mula for i? as a function of gs can be found in 
Eq. (7) of Ref. [10]. We notice that the depen- 
dence of R on gs is weaker than the dependence 
of B{B^ T^v) on gp, but the present la up- 
per bounds on the effective coupling constants 
are similar. However, if nature has opted for a 
large charged-Higgs contribution suppressing R 
to values below 0.2, it will not be easy to deter- 
mine gs because the curve in Fig. [3|is quite flat. 
At present there is also an open experimental is- 
sue in i? ^ DtVt'- The Monte Carlo simulations 



-fPor a lattice study of aY in the B ~> tt form factor, 
which involves the much larger range < |2| < 0.28 than 
our B ^ D transition, see Paul Mackenzie's talk at this 
conference. 
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Figure 3. i? EE B{B -> DTVr)IB{B -> Divf) 
as a function of gs- Light gray band: i?°''P = 
0.416 ± 0.117 ± 0.052 [21]. Dark gray band: R 
computed with BELLE form factor Vi [T7] . Gray 
band: R computed with HFAG form factor Vi 

m- 

use the SM formula for the decay distribution in 
B Drvr, but the D energy is higher on av- 
erage if a charged-Higgs contribution is present, 
because there is a destructive interference with 
the longitudinal W boson contribution. The effi- 
ciencies for D detection are smaller for very soft 
D's, which may affect the upper bound on gs de- 
rived from _R°'^P. 

It is well known that the sensitivity of i? ^ 
DtVt to charged-Higgs effects is improved, if in- 
formation on the T polarisation is included [55]. 
However, in B factories the r is too slow to decay 
with a displaced vertex, so that the r kinematics 
is not accessible to experiment. To my knowledge, 
the only theory paper addressing this problem is 
Ref. [55] , which proposes to study the differential 
decay rate dT/dEi), where Ed is the D meson en- 
ergy in the B rest frame. We have studied the de- 
cay chain B — > DDrT~ with the subsequent decay 
ti^Vt and propose to study the triple dif- 
ferential decay rate dT / {dED dE.^ d cos 6) , where 
Et^ is the 7r~ energy and 9 is the angle between 
the three-momenta of the D and the tt^. This 
quantity not only discriminates between SM and 
charged-Higgs effects in an excellent way, it also 
allows to constrain the phase of gs- This feature 
is illustrated in Fig. [3J 
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